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Abstract 

We calculate the cross section for the inclusive production of B mesons in pp and pp 
collisions at next-to-leading order in the general-mass variable-flavor-number scheme 
and show that a suitable choice of factorization scales leads to a smooth transition 
to the fixed-flavor-number scheme. Our numerical results are in good agreement 
with data from the Tevatron and LHC experiments at small and at large transverse 
momenta. 
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1 Introduction 


In the last twenty-five years, there has been much interest in the study of S-meson pro¬ 
duction in pp and pp collisions at hadron colliders, both experimentally and theoretically. 
First measurements were done in the late eighties by the UAl Collaboration at the CERN 
SppS collider |1] operating at a center-of-mass energy of VS = 0.63 TeV. Later measure¬ 
ments were performed by the CDF and DO collaborations at the Fermilab Tevatron in pp 
collisions at VS = 1.8 TeV piSlITIIB] and at y/S = 1.96 TeV piTlIS] . Recently, the CMS 
PIOIITT] . ATLAS pWS] . and LHCb [miT5] collaborations at the CERN Large Hadron 
Collider (LHC) published their results for inclusive and Bg meson production in 

pp collisions at y/S = 7 TeV. In all these measurements, the differential cross sections 
da/dpT, where pt is the transverse momentum of the B meson, integrated over specific 
rapidity {y) regions, or da/dy integrated over px > PT,min were obtained. While the CMS 
and ATLAS data were obtained at central rapidities, LHCb performed measurements at 
forward rapidities, 2 < y < 4 . 5 . Actually, only few measurements were able to explore the 
small-pT range down to PT,mm — 0, namely those by CDF [6] and the two by LHCb [ 141115 ] . 
In all other data, the cut PT,mm > 5 GeV was imposed. A unified theoretical description 
of these data, which covers both the very small and the large pt range, requires special 
efforts, which we shall undertake in this paper. 

On the theoretical side, it is generally accepted that for the treatment of R-meson pro¬ 
duction at small pr values, as well as for the calculation of the integrated cross section 
including the small-pr range, i.e., with pt of order 0{mb) and below, where mb is the 
5-quark mass, one should use the so-called massive scheme or fixed-fiavor-number scheme 
(FFNS) [T6lll7lll8l[T^ ?]. in which the number of active quark flavors in the initial state 
is limited to uj = 4, and the b quark appears with explicit mass dependence only in the 
final state. In this case, the b quark is treated as a heavy particle in the final state and 
not as a parton in the initial (anti)proton. 

In the large-pT region, characterized by pt ^ the so-called massless scheme or zero- 
mass variable-fiavor-number scheme (ZM-VFNS) [^22f23ll24|l25|l26ll27p28f29ll3n|l31|l32] is 
considered to be appropriate. This is the conventional parton model approach, where 
the b quark is considered massless like any other parton. The b quark is also treated 
as an incoming parton coming from the (anti)proton leading to additional contributions 
from hard-scattering subprocesses besides those with u, d, s, and c quarks or the gluon 
{g) in the initial state. Although this approach can be used as soon as the factorization 
scales Pi and pp associated with the initial- and final-state singularities are above the 
starting scale of the parton distribution functions (PDFs) and fragmentation functions 
(FFs), its predictions are reliable only in the large-pr region, pp 3> mi,, where terms of 
the order of m^/pp can be safely neglected. A next-to-leading-order (NLO) calculation in 
this scheme automatically resums leading and next-to-leading logarithmic (NLL) terms. 
At the same time, all non-logarithmic terms through 0{as) relative to the Born approx¬ 
imation are retained for mb = 0. With the conventional choice of renormalization and 
factorization scales, Pr = pi = Pf = mp with mp = Jpp + mf, the results are dominated 
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by contributions from the 6-quark PDF down to pt — 0. 


On the other hand, the general-mass variable-flavor-number scheme (GM-VFNS) 
provides a theoretical framework which combines the FFNS and the ZM-VFNS. The ZM- 
VFNS is extended into the intermediate-pr range by retaining the mass-dependent terms 
of the FFNS. With the conventional choice of scales indicated above, the results in this 
scheme are also dominated by the contributions of the 6-quark PDF down to pt — 0. As 
a consequence, there is no smooth transition from the GM-VFNS at large px values to 
the FFNS in the small-pr range, and the GM-VFNS fails to describe the small-pT GDF 
and LHGb data. 


The dominance of contributions with 6-quarks in the initial state at small px values is 
linked to the fact that this part is treated in the massless scheme, as a calculation of the 
6-quark-initiated subprocesses in a scheme with massive partons (like the AGOT scheme 
|47j ) is not available for hadropro duct ion The cross section with massless partons is, 
however, divergent for t 0. For a realistic description, we thus have to hnd a way 
to eliminate or modify this contribution in the small-pr region. In this paper, we shall 
develop an approach to modify the GM-VFNS in such a way that it matches the FFNS 
with the exact m;, dependence by a suitable choice of p/ and pp. We shall study how 
these modihcations can lead to a better agreement with presently available experimental 
data at small px values. 

The content of this paper is as follows. In Section 2, we introduce our strategy for the 
transition to the FFNS and compare our predictions with the GDF [6] and the LHGb [15] 
data. We shall also present results to be compared with measurements by the ATLAS 
Gollaboration [13]. Our conclusions are presented in Section 3. 


2 Small-px results and comparisons with data 

In this section, we shall discuss a viable unihed framework of theoretical predictions for 
inclusive H-meson production at small and large px values. We shall compare with the 
cross section distributions da/dpx measured by GDF |6| and LHGb [T5|. Throughout this 
paper, take the 6-quark pole mass to be rrib = 4.5 GeV, evaluate at NLO with 

Uf = A and A^ = 328 MeV if pp < mt, and with nj = 5 and A^ = 226 MeV if pp > mi,, 
and use the GTEQ6.6M proton PDFs [IH] unless otherwise stated. 

We start with results to be compared with the GDF data |6]. In Fig. we show NLO 
predictions in the FFNS with Uf = A. The full line shows the result for the default choice 
of scales, pi = ^imx with = 1 for i = R,I, while the dashed lines represent an estimate 
of the theoretical error obtained in the usual way, by varying up and down by a factor 
of 2. We take the transition of 6 and 6 quarks to the observed B mesons into account by 

^ For deep inelastic scattering, heavy-quark-initiated processes at NLO with massive quarks have been 
considered in Ref. [48] . 
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Figure 1: NLO predictions for da/dpr of pp —)■ + X with y/S = 1.96 TeV and 

—1 < 2/ < 1 in the FFNS are compared with CDF data |6]. The dashed lines represent 
the theoretical uncertainty estimated by varying pR and pi up and down by a factor of 2 
about the default choice Pr = pi = mR. 

using the branching fraction B{b ^ B) = 39.8% [50] as an overall normalization factor. 
The prediction in the FFNS agrees with the CDF data quite well, within experimental 
errors, up to pr — 15 GeV. Beyond this value of pr, the FFNS starts to overestimate the 
data, as has been shown already in our previous publication uni. 

In the FFNS, there is no need for FFs. However, a pir-independent FF might be introduced 
on phenomenological grounds and because of theoretical considerations to guarantee a 
proper matching between the schemes with Uf = A and n/ = 5. In the left panel of Fig. 
we show results obtained using the /ii?-independent Peterson FF [51] with parameter 
e = 10“^. We hnd only marginal differences with respect to the case where a constant 
branching fraction is used. Note that there are no g,q,q ^ B transitions in the FFNS. 

In addition to uncertainties from scale variations, there are also uncertainties due to errors 
in the input. We postpone the discussion of errors in the parametrizations of the PDFs to 
when we present predictions for the LHCb experiments, in Fig. below, but instead show 


4 





Figure 2: NLO predictions for da/dpT of pp — )■ + X with \/~S = 1.96 TeV and 

< ?/ < 1 in the FFNS are compared with CDF data |6]. Left panel: evaluation using 
the Peterson FF with e = 10“"^ and varying and pj by a factor of 2 up and down about 
the default choice Pr = pi = rriT- Right panel: uncertainties from variations of mi,. 

the influence of variations on the default prediction for the Tevatron measurements in 
the right panel of Fig. At small pt values, the uncertainty is comparable in size with 
the scale uncertainty, but it is negligible for px ~ 2mfe. 

In Ref. [10], we presented detailed comparisons of FFNS and GM-VFNS results with CDF 
data for production (see, e.g.. Figs. 7 and 8 in Ref. jlOl)- For calculations in the GM- 
VFNS, we use the pir-dependent FFs described in Ref. m. In this case, there are also 
small contributions due to transitions from light quarks and gluons to B mesons. With 
the default choice = ii = ip = 1, the GM-VFNS predictions diverge for pt —)■ 0, in 
obvious disagreement with the data. We notice, however, that the FFNS and GM-VFNS 
predictions approach each other at around pt — 20 GeV, i.e. 4 to 5 times mi,. In Fig. 
we show a similar comparison of the FFNS and GM-VFNS predictions for the case of 
the LHGb data [T5|. Although these data correspond to much higher y/S values and 
to different y ranges compared with the previous results in Ref. [lO], we observe similar 
qualitative behaviors of the FFNS and GM-VFNS predictions and a transition point at 
about the same value of pt, namely px — 20 GeV. 

One should expect that a correct treatment of kinematic constraints due to the hnite 
heavy-hadron mass is important, in particular at small px values. However, there is no 
unique prescription to take into account the hnite masses of the heavy quark and the 
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Figure 3: NLO predictions for da/dpT of pp —?■ B^+X with \/~S = 7 TeV and 2 < y < 4.5 
in the FFNS (dashed line) and the GM-VFNS (solid line) are compared with LHCb data 

m- 


heavy hadron at the same time. A prescription for including the heavy-hadron mass in 
connection with massless quarks based on the light-cone scaling variable was introduced 
in Ref. [52] • Here, we propose a generalization to the massive-quark case. This amounts 
to scaling the partonic cross section for the production of the massive quark a, with 
mass nia, energy Ea, and three-momentum pa, that initiates the formation of the heavy 
hadron H, with mass Mh, energy Eh = \lMjj + p^ cosh ?/, and three-momentum pn = 

Mjj sinh^ y + Pt cosh^ y, as 


da„ -)■ 


daa 

Rl 


Ra^H — 1 — 


Mfj - 


a^H 


{Eh+PH y - z'^ml ’ 


( 1 ) 


where z = {Eh +PH)/{Eh + Ph) is the light-cone scaling variable. In contrast to naive 
expectations, one hnds Rh < 1 corresponding to a slight enhancement of the cross section. 
In Fig. m we show the result of a calculation where this correction factor is taken into 
account. The effect is small mainly because it is only the difference between the b- 
quark and R-meson masses that enters. Therefore, also the additional suppression from 
tighter phase space limits is numerically not relevant. The poor small-pr behavior of 
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Figure 4: NLO predictions for da/dpr of pp —)■ + X with \/~S = 1.96 TeV and 

— 1 < 2/ < 1 in the GM-VFNS are compared with CDF data [6]. Besides the default 
prediction for scale parameters = 1 (solid line), a theoretical-error band encompassed 
by the predictions for = 1/2 (upper dashed line) and = 2 (lower dashed line) is 
shown. The effect of including quark and hadron mass corrections via the phase space in 
the default prediction is also shown, for = 5.28 GeV (dotted line). 


the GM-VFNS calculation cannot be remedied by such a naive treatment of phase space 
restrictions. 

The results shown in Figs. Bi suggest that the predictions of the FFNS and the GM- 
VFNS are very similar in the range 15 GeV ^ Pt ~ 20 GeV. In this range, both theories 
are justihed , with the FFNS naturally extending to smaller px values and the GM-VFNS 
to larger pt values. Therefore, one could designate a transition point pt at pr — 20 GeV, 
say, where one switches from one scheme to the other [53)154] . For practical purposes, this 
would require the knowledge of the PDFs in the scheme with n/ = 4 active flavors up to 
Pi = Pt [55]. A naive prescription to combine the two schemes could be to use matching 
functions 9{p^ — p\) and [1 — 6{p^ — p^)] to multiply the FFNS and GM-VFNS cross 
sections, respectively. Such a prescription will, however, lead to a discontinuity in the 
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Pt (GeV) 


Figure 5: Subtraction term for the gg channel evaluated with zero (solid line) and 
hnite (dashed line) value in the NLO prediction for da/dpT of pp —)■ + X with 

= 1.96 TeV and -1 < ?/ < 1 in the GM-VFNS. 


combined cross section at pt = Pt, a property that is certainly unphysical. One could try 
to invent different ways to combine the two schemes by introducing some smooth matching 
function, e.g., G{pt) = Pt/{Pt + with c = 5 [SH], and imposing, schematically, the 

prescription a = G(pt)o'gm-vfns + [1 ~ G'(pr)]o'FFNS- However, such an approach would 
introduce a new parameter, and it is unclear how theoretical uncertainties related to the 
choice of this parameter should be estimated. 

It will be helpful to take a closer look at the structure of the partonic cross sections in 
the GM-VFNS. They involve subtraction terms needed to separate the large logarithms 
ln(pf./m^) at hxed relative order For example, the cross section of the partonic subpro¬ 
cess gg —)■ gbb contains terms with large logarithms which can be written, schematically, 
as convolutions f{g —>■ 6) ® cr^gb —?• gb) and a{gg —>■ gg) ® d{g —)• b) with a partonic 
PDF f{g —)• b) and a partonic FF d{g —)■ b) (see Ref. [37] for dehnitions and a precise 
formulation of the corresponding cross sections). These splitting functions are absorbed 
in the PDFs and FFs, respectively. The large logarithms subtracted in this way are added 
back by calculating the cross sections of the gb —)■ gb and gg —)• gg subprocesses with a 
6-quark PDF and a g ^ B FF, respectively. There are similar subtraction terms and 






corresponding contributions involving PDFs and FFs of light quarks and the gluon in 
other channels as well. 

In Fig. we show numerical results for the subtraction terms needed in the gg channel 
(see Eqs. (45)-(53) in Ref. [37]), using a linear scale for better visibility. Their contribution 
is small, but not negligible, even at large pt values. Specihcally, the evaluations using 
the correct prescription with rrih 7 ^ 0 (dashed line) and the approximate one with mb = 0 
(full line) are compared with each other. We observe from this that such mass effects are 
small and cannot be responsible for the unphysical increase of the total result for da/dpT 
towards small pt values. 

It is, of course, unavoidable that the subtracted terms differ from those added back in 
the PDFs and FFs. The subtractions are obtained at hxed order, 0{as) in our case, 
while the PDFs and FFs contain the large logarithms resummed to all orders. Formally, 
the differences are of higher order in However, these higher-order terms are folded 
with cross sections calculated in the ZM-VFNS, and these cross sections are singular for 
Pt —t 0. Therefore, it is not surprising that the contributions with b quarks in the initial 
state dominate at small pr values. 

Obviously, NLO cross sections of 6 -quark-initiated subprocesses that are convoluted with 
PDFs and FFs evolved at NLO contain terms singular for —)• 0 at one order beyond 
the subtracted terms. In fact, it would be a major task to derive the missing next-to- 
next-to-leading-order (NNLO) subtraction terms. With such an extended version of the 
GM-VFNS including NNLO subtractions, but still using partonic cross sections derived 
in the ZM-VFNS, the problem would be shifted to one order higher, but it remains to 
be seen whether the required cancellations of singular terms can be obtained with the 
required numerical precision. In addition, strictly speaking, the NNLO subtraction terms 
only make sense in combination with the hxed-order calculation at NNLO. Otherwise, at 
large pt values, the NNLO subtraction terms and their NLO hxed-order counterparts do 
not cancel, and the ZM-VFNS is not recovered. 

In turn, one could argue that 6 -quark-initiated processes evaluated with 6 -quark PDFs 
and FFs should be included in the GM-VFNS at LO only. The predictions thus obtained 
EH exhibit better agreement with data in the medium-pr range, between 2 and 7 GeV 
or so, but switching off these NLO terms is again not sufficient to completely eliminate 
a singular behavior for pt —t 0. Moreover, and most importantly, numerical evaluations 
show that NLO corrections in the zero-mass part of the GM-VFNS are essential at large 
Pt values. Therefore, we do not follow this option either. 

Instead, we try to exploit the freedom offered by the presence of pr, p/, and pf, pa¬ 
rameters that are present anyway. Their values are not determined by theory, but some 
choice has to be made, based on some reasonable, but ad-hoc physical argument. In fact, 
a judicious choice of scales can lead to a suppression of the potentially dangerous con¬ 
tributions from initial-state 6 quarks. This exploits the fact that all commonly available 
PDF hts assume that the 6 -quark PDF is zero below some starting scale, usually chosen 
to be Pi = mb. The same is true for the FFs: the FF for the b ^ B transition vanishes 
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Figure 6: NLO predictions for da/dpT of pp —)■ + X with y/S = 1.96 TeV and 

—1 < 1/ < 1 in the GM-VFNS are compared with CDF data [B]. Left panel: evaluations 
with = 1 and = 0.5, (solid line), 0.4 (lower dashed line), and 0.6 (upper 

dashed line). Right panel: evaluations with = 0.5 and = 1 (solid line), 0.5 

(upper dashed line), and 2 (lower dashed line). If ^i^f < 1, then pj^F = '^b is put for 
Pt < - 1. 


for /i^ < mb. Therefore, with pi p = a value ^j^f < 1 will render the 6-quark 

PDF and FF zero for pt < mb^l/i'j^p — 1. 

In Fig. 1^ we show results obtained with such scale choices and compare them with CDF 
data |6]. Indeed, values for and ^f of about 1/2 lead to the required suppression of 
6 -quark-initiated contributions. Specihcally, in the left panel of Fig. we choose = 1 
and = ^F = 0.4,0.5, 0.6. In all these cases, there is a turn-over near pt = 2.5 GeV, 
and the agreement with the CDF data is reasonably good. We impose the additional 
constraint that pi and pf are not allowed to take values pi,pF < nT^b, i-e. the DGLAP 
evolutions of the PDFs and FFs are frozen below this scalej^ This explains the slight 
bumps that occur in the pt distribution at pt = 6.0, 7.8, and 10.3 GeV for ii = ^f = 0.6, 
0.5, and 0.4, respectively. Obviously, the freedom in the choice of the default values of 

^In fact, we have to freeze /i/ slightly below mb, pi = Cmb with C < 1, since the 6-quark PDF 
parametrization of the CTEQ6.6M set, which we use here, vanishes only strictly below threshold. For 
our numerical evaluations, we choose C = 0.99. 
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^7 and leads to a moderate extra uncertainty close to the transition region, but the 
results in the small-p^ range, pr ~ 5 GeV, and for large px values stay unaffected. 

Taking now = 1, = ^f = 0.5 as the new default scale choice, we estimate the 

theoretical errors in the usual way by varying the scale parameters by a factor 2 up 
and down about the default scale choice. In fact, it turns out that the variation of pr 
is the dominant source of the theoretical uncertainties, and we simplify the subsequent 
calculations by only considering variations of ^r in the range from 0.5 to 2. Note that we 
do not introduce an extra prescription to freeze pr below mb because, hrst, the choice of 
Pr is not related to switching off 6 -quark-initiated subprocesses and, second, full variations 
of pr are needed to obtain realistic estimates of the theoretical uncertainty. The resulting 
error band is shown in the right panel of Fig. We emphasize that the freedom in 
choosing specihc values of ^7 and ^f as default does not introduce a large additional 
uncertainty, as may be understood by comparing the left and right panels of Fig. 

In Fig. we compare the LHCb data [15] with our theoretical predictions implemented 
with the scale conventions introduced above. In the left panel, we see again nice agreement 
between the FFNS predictions and experimental data. The right panel tells us that the 
data are not quite as well described by the GM-VFNS predictions, but the agreement is 
quite acceptable, as the data points are covered by the error band for px ~ 2.5 GeV. In 
this px range, the agreement of the data with predictions obtained in the FONLL scheme 
Ba is quite similar na. 

In Fig. we show a comparison of results obtained using different PDF parametrizations. 
Specihcally, we use the GTIO [58], MSTW [59], HERA-PDF1.5(NLO) [60], and NNPDF 
2.3 [6T1I62] PDF sets. The GTEQ 6 . 6 M [19] PDFs adopted elsewhere in this paper yield 
a result that is very similar to the one obtained using the GTIO PDFs. We observe that 
there are differences in the small-p^ range, px ^4: GeV, which are somewhat larger than 
the experimental errors. We should, therefore, expect that R-meson production data at 
the LHG will help us to further constrain the PDFs. In particular, there is sensitivity to 
the gluon PDF in this kinematic range. 

For completeness, we also consider the production of and mesons. Appropriate 
experimental data were published by the LHGb Gollaboration in Ref. [15]. In Fig. we 
present comparisons with NLO predictions in the GM-VFNS using the scale setting and 
theoretical-error estimation prescriptions described above to hnd good agreement. 

Finally, we compare i?’''-meson production data taken by the ATLAS Gollaboration [13] 
with our NLO GM-VFNS predictions. These data extend into the very-large-pr range, 
where we expect the GM-VFNS to be quite appropriate. Indeed, we hnd good agreement, 
except for the lowest px bin, 9-13 GeV, and for central rapidities, where the data are 
slightly overestimated. Because of the large px values probed, the adjustment of scales to 
match to the FFNS as described above is not an issue here. 
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Figure 7: NLO predictions for da/dpT of pp —?■ B^+X with \/~S = 7 TeV and 2 < y < 4.5 
in the FFNS (left panel) and the GM-VFNS (right panel) are compared with LHCb data 
[T5] . The default scale choice is = 1 and 0 = 0.5 (solid lines), and the theoretical- 

error bands are obtained by varying by a factor of 2 up (lower dashed lines) and down 
(upper dashed lines). If < 1, then = mi, is put for pp < — 1. 


3 Conclusions 


Any theoretical prediction for hadronic collisions within perturbative QCD requires the 
factorization of initial- and hnal-state singularities. This unavoidably introduces factor¬ 
ization scale parameters, which cannot be predicted from theory. Their choice must be 
based on physical arguments. We exploited the freedom in this choice to hnd a pre¬ 
scription that extends the reliability of NLO predictions in the GM-VFNS down to small 
Pt values. With scale parameters .^/ < 1, it is possible to eliminate contributions from 
the heavy quark in the initial state. These contributions, dominated by the subprocess 
gb —)■ gb^ are treated in the parton model with zero quark masses and would lead to an 
unphysical increase of the cross section for t 0. We showed, however, that with a 
judicious choice of and one can switch off these contributions in the small-p^ 

range. Our prescription leads to a modihed GM-VFNS yielding results in good agreement 
with recent data from the Tevatron and LHG experiments. 
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Figure 8: NLO predictions for da/dpT of pp —?• with y/S = 7 TeV and 2 < y < 4.5 

in the GM-VFNS with default scale choice are compared with LHCb data [I5]. The PDF 
error is estimated by using the CTIO (solid line) [58], MSTW (dashed line) [59], HERA- 
PDF1.5(NLO) (lower dotted line) [60], and NNPDF 2.3 (upper dotted line) [6T1I62] sets. 
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Figure 9: NLO predictions for da/dpT of pp — B^/B^ + X (left panel) and pp —?• 
Bg/B^g + X (right panel) with \/S = 7 TeV and 2 < y < 4.5 in the GM-VFNS are 
compared with LHCb data ng. The default predictions and theoretical-error bands are 
evaluated as in the right panel of Fig. 
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Figure 10: NLO predictions for da/dpT of pp —)■ B~^ + X with \/S = 7 TeV and 
0 < 7/ < 0.5 (upper left panel), 0.5 < y < I (upper right panel), 1 < t/ < 1.5 (lower 
left panel), and 1.5 < y < 2.25 (lower right panel) in the GM-VFNS are compared with 
ATLAS data ra. The default predictions and theoretical-error bands are evaluated as in 
the right panel of Fig. 
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